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1. Introduction 



About sixty years ago, Newton and Wigner showed that there is no position operator 
for a massless particle with spin higher than 1/2 [1]. This fact lead many authors to 
conclude that it is not possible to define a wave function for a photon, which has zero 
mass and spin (or helicity) 1. However, Bialynicki-Birula |2][3] and, independently, Sipe 
[I] , introduced a function of the position and time coordinates that completely describes 
the quantum state of a photon. In their opinion, which is shared by us, such function 
may be referred to as the photon wave function. The wave equation for this function 
can be derived from the Einstein kinematics for a particle with spin 1 and zero mass 
in the same way that the Dirac equation is obtained for a massive particle with spin 
1/2 [51 El- A strong argument in favor of this photon wave function formulation is that 
the corresponding wave equation is completely equivalent to the Maxwell equations in 
vacuum. 

In the present work we extend the Bialynicki-Birula-Sipe formalism to include 
the interaction of photons with non-absorptive continuous media. In section [2] we 
briefly discuss the photon wave equation in vacuum and some properties of photon 
wave functions. In section |3] the photon wave equation is modified including a term 
associated with the induced current density in the medium, that makes it equivalent to 
the macroscopic Maxwell equations for electromagnetic waves in the presence of matter. 
In section 0] we apply the second quantization procedure to the photon wave function 
formalism in a linear non-dispersive and non-absorptive medium, allowing for description 
of quantum phenomena. In section [5] we present a method for treating the scattering of 
photons by material media. As an example of application of the proposed formalism, 
in section |6] the state of the photon pairs generated by parametric down conversion is 
obtained with an intuitive derivation. In section [7] we present our concluding remarks 
and summarize our results. 



2. Photon wave equation in vacuum 



The Bialynicki-Birula-Sipe wave equation for a photon described by the wave function 
*(r,t) is |2l Eli 

ih— — = hcaV x (1) 
dt w 

where % is the Planck's constant divided by 2n and a is an operator that acts on the 

photon helicity. A spin one can be represented by a vector in the 3-dimensional complex 

space, so the vectorial character of \P incorporates its internal degrees of freedom. The 

photon wave function can be decomposed in the eigenvectors of the helicity operator a 

and also in its real and imaginary parts: 




* = * + + with * ± = J— E± ± h B ± and a* ± = ±* ± . (2) 
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If we associate the real vectors E = E + + E_ and B = B + + B with the electric and 
magnetic fields of the photon and e$ and /xq with the permittivity and the permeability of 
free space, (CEJ) together with the divergence-free condition of the photon wave function 
V • Vl/ = (which is essential to arrive at the photon wave equation from Einstein 
kinematics E]) are completely equivalent to the set of the four Maxwell equations 
in vacuum. The eigenvalues of the helicity operator correspond to photons with left 
and right circular polarizations. An arbitrary polarization state can be written as a 
particular combination of these two. 

There is a peculiarity in the form ([2]) for the photon wave function: Its modulus 
squared has the unit of energy density, not of probability density, and agrees with the 
classical electromagnetic energy density. This fact is associated with the problem of 
non-localization of the photon (Tj, and it happens that the energy density is the local 
useful concept, since photons have no other scalar quantities such as mass or charge, but 
only energy [1]. The normalization condition is that the integral of \^\ 2 is equal to the 
expectation value of the photon energy. There is no direct association of the modulus 
squared of the wave function with the probability density of detecting the photon in the 
region, unless we consider dimensions larger than the photon wavelength [7] . But neither 
the localization of massive relativistic particles is a clear concept. In this sense, photons 
are not significantly different from massive particles, as both can have an exponential 
falloff of their wave functions with distance in relativistic treatments [Bj, the difference 
being that photons are always relativistic. 



3. Photon wave equation in material media 

Our proposal to treat the interaction of photons with continuous media is to include a 
term proportional to the current density J induced in the media due to the presence of 
the photon in (JTJ: 

9* J 

ih— = UcaV x * - ih—= . (3) 

at v 2eq 

This approach is similar to the one adopted by Keller j9]. However, Keller was interested 
in calculating the wave function of a photon emitted from an atom and the current 
density was a source term. Here we are considering continuous media and the current 
density J represents a response of the medium due to the presence of the photon, not 
being an independent source of photons. So the objective of the present work is very 
different from Keller's, although the initial step is similar. When considering continuous 
media, we are implicitly assuming that the photon wavelength is much larger than the 
distances between the atoms that compose the medium, in such a way that the collective 
motion of the charges may be approximated by a continuous current density. We can 
write 

J = J f + VxM + dP/dt, (4) 

where M and P are the magnetization and electric polarization of the medium and Jf 
is the free current density, all of them being functions of the electric and magnetic fields 
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of the photon. For instance, in a non-dispersive dielectric medium the polarization can 
be written as 

P,(r,t) = 5>g } J5,-(r,f) +'£j$ k E J (r,t)E k (T,t) + (5) 

3 3,k 

expanded in linear and non-linear responses to the photon field. Since the non-linear 
response terms represent situations in which there are more than one photon involved, 
a procedure of second quantization of the photon wave function is necessary to correctly 
describe them. This issue will become clear in section [6] where we treat the parametric 
down-conversion, which arises from the second-order non-linear response of the medium. 
So, if we want to deal with the wave function of a single photon, we must disregard the 
non-linear response of the medium. In section [6] we will be interested in the interaction 
of photons with non-linear media after applying the second quantization procedure to 
the photon wave function, so we will continue referring to the current density J in this 
section without specification of its linear or non-linear behaviors. 

Using the photon wave function from (J2J), we can see that (j3J is equivalent to 
the Maxwell equations by writing its real and imaginary parts. Its real part with the 
application of the operator a and its imaginary part give: 

VxE_, (6) 

B <9E <9D 

V x — = £ o^T + J , or V x H = J f + — , 7 

fjio at ot 

with E = E+ + E , B = B+ + B , H = B//i - M and D = £ E + P. Taking the 
divergence of the above equations we obtain 

|(V-B) = 0. (8) 

_(V-D) = -V-J f , or-(V-D- Pf ) = 0, (9) 

where pt is the free charge density and we have V ■ Jf = —dpi/dt. What the above 
equations tell us is that if V ■ B = and V ■ D = pf at any particular instant of time, 
which is true before the photon interacts with the media, these relations must be always 
satisfied. With this argument we complete the proof that (J3]) is equivalent to the set 
of Maxwell equations for electromagnetic waves. So our treatment predicts that the 
behavior of photons in linear non-absorptive media, where a one-photon theory may be 
used, is completely analogous to the behavior of classical electromagnetic pulses. This 
statement is supported by experiments [10] and is difficult to obtain with the formalism 
of second quantization JTTJ [121 131 EEU EES]. So we believe that in the same way that 
the Schrodinger equation is used to treat the behavior of massive particles under the 
action of a classical potential, the formalism of this section may be used to describe 
the interaction of photons with material media when there is no absorption or emission 
of photons. In cases where there is absorption or emission of photons we must apply 
the second quantization procedure on the photon wave function, to be done in the next 
section. 
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From ([3]) we can derive a continuity equation for the photon wave function. Taking 
the scalar product with on both sides of ([3]) and summing the complex conjugate, 
after verifying that (crV x \I>)* • \I> = (V x • and using the vectorial identity 

V • (A x B) = B • (V x A) - A • (V x B), we arrive at 



where S = — ic^* x (<r\l/) is the current density of the photon wave function. Using 
(J2]) and taking the time average of S over one period of oscillation of the photon field, 
we have (S) = (E x B)//xq, analogous to the electromagnetic energy flux. The term 
(\|> _|_ vEr*) . J/y/2eo = E ■ J represents the density of work done on the charges of the 
medium per unit time. The continuity equation for the photon wave function is an 
energy continuity equation. In our formalism, we are assuming that the photon energy 
density has always the same form when written in terms of the photon electric and 
magnetic fields, independently of the medium, and that the material energy density is 
calculated by the work done on the charges of the medium, without distinction whether 
the work is done on bound or free currents. This approach differs from the one adopted 
by Bialynicki-Birula to treat photons in material media j3], where the photon wave 
function is written as F = D/ y2e ± iB/ y/2jl and there is no current term in the photon 
wave equation, since the bound currents are included in the photon wave function via 
D, Xe and Xm- The photon energy density |F| 2 is a function of the electric permittivity 
e and magnetic permeability fi of the medium since we have D = eE and H = B//i in 
a linear isotropic medium. 

When a photon enters a medium, part of its energy is transferred to the medium. 
The same occurs with its momentum. The discussion of the correct form for the 
momentum and energy densities of classical electromagnetic waves in a medium has 
been the subject of a long debate |16j . The eventual conclusion is that there are 
many compatible ways of defining the electromagnetic and material parts of the energy 
and momentum densities of electromagnetic waves in material media. Here we are 
assuming a particular division in the quantum domain, where the electromagnetic parts 
are |^| 2 = £ |E| 2 /2 + |B| 2 /(2/i ) for the energy density, S for the energy flux and S/c 2 
for the momentum density. This division has been successfully tested with reasonable 
models in the classical domain [TTt [T8] . 

Our formulation gives the same results as Bialynicki-Birula's for linear media: 
photons obey the Maxwell equations, the difference being that we implicitly divided 
the total energy and momentum of the photon in the medium in their electromagnetic 
and material parts, associating only the electromagnetic part with the photon wave 
function. As a consequence, our formulation can be readily extended to non-linear media 
after the procedure of second quantization of the photon wave function, representing an 
advantage to treat more general cases, as we discuss in the next sections. 




VS 




1 



(10) 
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4. Second quantization of the photon wave function in a linear 
non-dispersive and non-absorptive medium 

The second quantization procedure can be applied to the one-photon theory to build a 
quantum field theory for describing states with an arbitrary number of photons. This 
procedure is described in [6]. Here we will make a similar second quantization procedure 
in an infinite linear non- dispersive and non-absorptive medium. In such a medium, the 
material acquires an energy density u mat = Xe^o|E| 2 /4 + x m |B| 2 /[4(l + x m )/io] due to the 
presence of the photon [18] , Xe and Xm being the electric and magnetic susceptibilities of 
the medium, apart from the electromagnetic energy density \^f\ 2 = £o|E| 2 /2+|B| 2 /(2/zo). 
So the total energy density of the photon in such a medium can be written as |\E ,/ | 2 , 
where = *Sf' + + *Sf'~ is the "dressed photon" wave function, with 

* 7± = y^E± ± i^B ± and a*'* = ±V ± , (11) 

where e = (1 + Xc)^o an d A* = (1 + Xm)/-*o- It is straightforward to show that 
| | 2 = |v|>| 2 -)- u mat . This form for SI/' is equivalent to the one used by Bialynicki- 
Birula in the case of non-dispersive linear media [3]. However, if dispersion is included 
the modulus squared of the Bialynicki-Birula's form for the photon wave function in the 
medium does not give the total energy density anymore, while the division of the total 
energy density in electromagnetic and material parts remains valid [TS]. Although we 
will not treat dispersive media here, it is important to stress this fundamental difference 
between the treatments. The above relations can be inverted, giving 

E = E+ + E_ = -= (*' + *'*) , (12) 
V2e 



B = B + + B = -i^f [(a*') - (a*')*} ■ (13) 

The total energy of the "dressed photon" can be written as 

H = J d 3 r|*'(r,t)| 2 . (14) 

The "dressed photon" wave function can be decomposed in its energy eigenstates, which 
are plane- wave modes: 



*'(r, t) = £ / d 3 ^^ a ks e^— ') e ks (15) 
with J2 S I \dks\ 2 — 1; where the plane- wave modes are defined by their wavevector k 



and polarization vector e ks , with frequency given by uo = ck/n, n = y (1 + % c )(l + Xm) 
being the refraction index of the medium. Substituting (fT5| in f lT4|) . the energy can be 
written as 

H = J2j <^khu\a^ s \ 2 . (16) 
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The above Hamiltonian is equivalent to the Hamiltonian of a continuous set of 
harmonic oscillators [19], and the result of the second quantization procedure is to raise 
the "dressed photon" wave function ( fl5|) to the status of a field operator 



*'(r, t) = £ / d 3 ^^ a ks e ks , (17) 

where ay_ s is the annihilation operator for photons in the mode (k, s) that, together with 
the creation operators a[_ s , obey the commutation relations 

[a ks , aw] = [aL» %] = °> [«k s , a w ] = <5 3 (k - k')5 aa /. (18) 
According to f[T2"l) and ([TBI , the electric and magnetic fields also become operators' 



E(r,t) = 2 / d 3 *J-?|-p e'< k "-"" a k ,e kJ + h.c, (19) 



2e(27r 



2(2^)3 - - [I X §ks J + (20) 

Equations ffT9|) and fl20|) for the electric and magnetic field operators are identical 
to the ones obtained with the traditional second quantization of the electromagnetic 
field in vacuum if we substitute e by e and \i by /x [H] . In the treatment of Smith and 
Raymer for the quantization of the photon wave function in vacuum [6], the operator ^ 
they find is equivalent to our 4f' + where 4?' is given by (fT7|) with Xc = Xm = 0. We 
believe that our form, in which has only annihilation operators, is more convenient, 
since it agrees with the condition 

* = & + »& < 21 > 

derived from ( II ip . 

A one-photon state can be written as 

|^W) = E/d 3 ^f)(k)aL|vac), (22) 

s J 

with J2 S J d 3 k\^'s s )] 2 = 1 and |vac) representing the vacuum state of the 
electromagnetic field. The corresponding one-photon wave function is [B] 



*' {1) (r,t) = (vacl^Cr.t)!^ 1 )) = E/ d3k \j^ys V'PW^-^s, (23) 



in agreement with (fl5l) . 

A two-photon state can be written as 



|^(2) } = E / d 3 fc | d'k' *J 2 J(k, kOaLaLy |vac> (24) 

s,s' 

X To arrive at (|20p from (fT5)) . we must use the fact that cre^ — i(k/fc) xe^, where is a polarization 
vector with helicity a. This can be verified choosing k/fc = z, et+ = (x+zy) / y/2 and §k- = (x— iy) j y/2. 
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and the corresponding two-photon wave function is 0, [2U] 

* ,(2) ( ri ,r 2 ,t) = (v a c\&( ri ,t)*'(r 2 ,t)\¥V) 

d 3 k r d 3 k' 



_ (2vr)3/2y (2tt)3/ 2 
x Vftwe^-^ek, (8) v^e^'^-^ew. (25) 

The tensor product in the above equation indicates the composition of two one-photon 
space states. Note that the wave function is automatically symmetrized via the term 
ty s J, (k, k') + \P S , s (k', k) . The extension to higher dimensions is straightforward. 

5. Scattering of photons 

Here we want to find the behavior of the operators E and B from ( II 9p and (1201) 
when, besides the linear medium with electric and magnetic susceptibilities Xe and Xm 
occupying all space, there is another medium in a finite region with different properties. 
Our treatment can be seen as an approximation of Keller's treatment for the photon 
scattering [9] [21] in the far-field regime. Applying the second quantization procedure 
including the interaction with matter, the operators ^ will obey the same wave equation 
([3]) as the wave functions So the operators E and B will also obey ([6]) and (J7J). 
Equation (CO) can be written in the following way 

Vx - = ^ + J'' ( 26 ) 
jx at 



where 



J' = J — V x 



XmB 



(1 + Xm)^o. 



dt 

contain terms due to the difference of the behavior of the scattering medium in relation 
to the quantization medium, according to ([!]). Because the current density J' is a 
function of the electric and magnetic fields, when the fields become operators J' also 
becomes an operator. 

Differentiating (125]) with respect to time, using ([6j) and the vectorial identity 
V x (V x E) = V(V ■ E) — V 2 E, we obtain a wave equation for the operator E: 

„2 A2"p 3]/ 

V 2 E - - — = -4vrP , with f c = -J*— - V(V • E). (28) 

The above equation is a wave equation in which f e is the source. An analog equation 
can be found for the operator B in a similar way: 

V 2 B - = > with f b = -V x J'. (29) 

c 2 at 2 /i 

Solutions of (|28p can be obtained using Green functions [22] . If we want to study 
the scattering of a photon in a known state in the paraxial regime, where the dependence 
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of the polarization with the propagation direction can be disregarded and a scalar theory 
may be used, the relevant solutions of fl28|) for the operator E(r, t) are 

E out (r,i) = E in (r,t) + / dV f dfG< + >(r, t; r', t')f c (r', t% (30) 

with G (±) (r,t;r',t') = — - / dA;^ -e^*"*), (31) 

on 2 j-oo |r — r'| 

where u = kc/n. 

Integrating fl3T|) in k, we see that the Green functions can be written as 

G (±) (M; r',0 = ° *(|r - r'| =F (£ — f)c/n). (32) 
47r|r — r | 

It can also be shown, by a direct integration on the variables <ft an d 9 in spherical 
coordinates ((f), 9, k) [E], that 

^fc e ip.(r-0-(t-f)] =G ( + )( r ,t ;r ' ; t -G(-)(r,t;r',t'). (33) 

If we are studying the scattering of a photon, the observation time t will always be 

greater than the time of the interaction t' . So, according to ( )32|) . (r, t; r', t') = 
and we can write 

GW(r,i;r',0 = / &W«0-«M)]. (34) 

(27T)' 3 7 fc 

The utilization of the Green function method to predict the behavior of the electric 
field operator, according to f )30|) . is very useful for calculating the state of scattered 
photons by a localized material medium. In linear media, we can use the same equation 
substituting the operator by the electric field of the photons, obtaining equivalent results 
to the scattering of an electromagnetic wave, since (|30|) is obtained from the Maxwell 
equations. But the interesting results come from the treatment of problems that cannot 
be obtained from classical electrodynamics. In the next section, we will apply the 
treatment to the scattering of a photon by a non-linear medium that can generate a 
photon pair in an entangled state. We will show that the quantum state of the generated 
photon pair we obtain is the same as the one previously obtained [23l HHJ [24] , but with 
a more intuitive interpretation. 



6. Treatment of parametric down conversion 

The parametric down conversion corresponds to the conversion of one incident photon in 
two by a non-linear birefringent crystal. As an example of application of the formalism 
developed in the previous section, we will calculate the quantum state of the photon pair 
generated in the process of type-I parametric down conversion, in which both generated 
photons have polarization orthogonal to the incident photon. To simplify the treatment, 
we will consider that the non-linear crystal is immersed in a medium with the same linear 
responses as the crystal. So the medium where the field is quantized occupies all space, 
and the non-linearities of the crystal will be responsible for the scattering processes. 
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The magnetic responses of these crystals are very small, so we can consider a refraction 



index n = J(l + Xe)- We will also disregard the dependence of the refraction index 
with the propagation direction and the effects of dispersion. These approximations are 
consistent with the observation of photons at small angles in relation to the propagation 
direction of the incident photon and to the use of interference filters on the detectors, 
post-selecting states with relatively narrow frequency spectra. 

With all these considerations, the non-linear components of the polarization of the 
crystal can be written as 

Pr(r,t)=Y,X%E J (r,t)E k (r,t), (35) 

where terms beyond the second order were disregarded. In type-I parametric down 
conversion, the generated photons have ordinary polarization and the incident photon 
has extraordinary polarization. According to (j4j), the non-linear current density 
associated with the ordinary polarization can be written as J' a NL = dP^ L /dt, and 
the scattering term, according to ( 1281) . is fl = —fXodJ' NL /dt. So the relevant source 
term for the scattering in (1301) is 

pp. pNL f)2 

r o (r>, t>) = -fM^- = -IMt-jjps [x£Kr')E (r>, t')E e (r', t')} , (36) 

where the subscript e (o) denotes extraordinary (ordinary) polarization and we have 
Xooe = Xooe + Xoel) si nce E (r',t') and E e (r',t') are operators that act on modes with 
orthogonal polarization and therefore commute. According to ( |25|) . if we denote the state 
of the incident photon by )> the two-photon component with ordinary polarization 
of the output state is 

V%(r u r 2 ,t) = (™\*Uo(ri,t)*Uo(r2,t)\1>®) . (37) 

Considering a first order scattering, in which one of the operators ^f' outo {r,t) above is 
substituted via (121]) by the scattered electric field of fl30|) and the other is kept with the 
non-scattered field, we can write 

^( ri ,r 2 ,t) = (vac\&M,t)n^E on ^,tM {i) ) (38) 

plus a term with the permutation ■<->■ r 2 , that correctly symmetrizes the wave function. 
From now on we will omit this permutation term to simplify the notation. The refraction 



index n appears in the above equation because we have yfe = y (1 + Xc)^o = n^/eo. 
Eonto is given by (1301) . with f e given by ( 1361) . 

The state \ip^ 1 ') has one photon with extraordinary polarization and in f )38|) 2 
operators that act on ordinary polarization modes are applied to the state |^ )■ We 
can substitute these operators by their commutator plus terms that result in zero when 
applied on \ip^). Using ffTTj) . (fl8|) and ( TT9|) . the relevant commutator to the problem 
in the limit that we disregard the dependence of the polarization with the propagation 
direction is [19] 



4f'( T , ti) J—E (V t') = [ d 3 k e i[k-(r-r')-c(t-tQ] 

W ° lri ' tlj 'V 2 ol ' j J J 2n(2vr)3 e 
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-G( + )( ri ,t;r',t'), (39) 



4c 2 dtdt' 

where (134")) was used in the last equality. 

Inserting the value of ( 1361) for f e in E outo (v2,t), obtained from (1301) . in ( 1381) and 
using fl39|) to substitute the product of the operators ^(r^ti) and Jeq/2E (t' ', t') by 
their commutator plus terms that result in zero when applied on \ip^), we have 



E' 3 (r\t')\ , (40) 



X dt' 2 



rj2 



with E' 3 (r',t') = (vac|nE e (r',t / )|V' (<) >. Substituting (El) and defining ^(r',*') 
I (i 3 fc 3 i?3(k3) exp(ik 3 • r' — zu^i), we have 

*foVi,r 2 ,t) = #^ [ dV f dt' / d% f d% [ d%^-^^l x 



(2n) e J ■/_«, 7 V J h k 2 
x ^'(^(r^^^kg^e 1 ^ 3- ^ - ^)' 1 ' -1 ^ 3- " 1- " 2 ^^ 1 ^ 1 ^ 1- " 1 *^ 1 ^ 2 ^ 1- " 2 ^ . (41) 

If we observe the generated photons at a distance from the crystal such that the time 
light takes to go from the crystal to the detectors is greater than the duration of the 
incident pulse, the integral in t' may be extended until t' = oo, resulting in a term 
proportional to 5(u3 — u)\ — c^)- So the term (0J3 — uj\) 2 in the above equation may 
be substituted by <J\. The above form for the wave function of the twin photons 
generated by parametric down conversion is equivalent to previous treatments that 
consider a Hamiltonian evolution of the quantum state of light interacting with a non- 
linear medium using perturbation theory [231 IT9] . 

If we use the form from (l3~TT) to the Green function, we arrive at a more intuitive 
form for the quantum state of the twin photons: 

*i 2) (ri,r 2 ,t) = - / dV/^dt'/ &k x J dk 2X ^ e (r')E' 3 (r',t') x 

x t { "f— { n • ( 42 ) 

|ri — r I |r2 — r I 

We can interpret the above wave function as being the result of the coherent 
scattering of one photon in two by all non-linear scatterers of the crystal. In each 
one of these processes, the incident photon generates a pair of photons in spherical wave 
modes emitted from each non-linear scatterer. The energy and momentum conservation 
in the process may be seen as consequences of the fundamental indistinguishability of all 
possible instants of time and positions in which a photon pair could be generated, which 
forces us to coherently add all the probability amplitudes. These probability amplitudes 
interfere destructively when energy is not conserved. However, although the sum of the 
energies of the twin photons is the energy of the incident photon, the energy of each 
photon can assume a wide range of values. So, the photon pair is highly entangled in 
energy. The probability amplitudes also interfere making the sum of the momenta of 
the twin photons be associated with the momentum of the incident photon, depending 
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on the crystal geometry since the crystal may absorb part of the momentum in the 
process. The momentum of each photon can also assume a wide range of values, so the 
generated photons are also entangled in linear momentum. 

Comparing ( 14"T|) and ( 125|) . we see that the post-selected wave function of the twin 
photons on the wavevector space is proportional to 

^(k 1; k 2 ) oc ^ / df / dV / d 3 fex£(r') x 

UiH2 J J J 

x ^ 3 (k 3 )e i ( k 3- k i- k 2)- r '-H"3-^ 1 -^)t'_ (43) 

Let us consider that the incident photon is in a beam mode that propagates in the z 
direction and that the crystal is a rectangular box with dimensions much larger than the 
beam waist in the x and y directions and a small width in the z direction. Let us also 
treat the collinear regime, in which the crystal is oriented to maximize the production 
of photons with wavevectors near the z direction. So the integrals in x and y in (1431) 
may be extended to infinity and the integral in z results in a constant. To illustrate 
the energy and momentum entanglement of the state, let us represent the vectors k 
as k = q + k z z and E 3 (k 3 ) as E 3 (u3, q%), since U3 and q3 determine k 3 . Let us also 
substitute the integral in the z component of k 3 by an integral in U3. With all these 
approximations, the integral in t' gives 6(003 — u\ — u 2 ) and the integrals on x' and y' 
give S 2 (q 3 — qi — q2). So the 2-photon wave function on wavevector space can be written 
within the paraxial domain (|q| <C |k|) as [21] 

*™ Qi) w 2, Q2) oc ^/ujiuj 2 E 3 (ut + u 2 , qi + q 2 ). (44) 

We cannot write this wave function as a product of a wave function for photon 1 and 
a wave function for photon 2, what characterizes the entanglement of the state in both 
energy and momentum. 

7. Conclusions 

We extended the Bialynicki-Birula-Sipe photon wave function formalism to include the 
interaction between photons and continuous non-absorptive media. The photon wave 
equation becomes equivalent to the macroscopic Maxwell equations in material media. 
Applying the second quantization procedure to the photon wave function, we obtained 
an alternative treatment for the quantum interactions between light and continuous 
media, that can be specially useful to treat the scattering of light by matter. As an 
example of application of the proposed formalism, the quantum state of the twin photons 
generated by parametric down conversion was obtained in agreement with previous 
treatments, but with a more intuitive interpretation. 
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